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ABSTRACT 

We present a calculation of the total cross section for top quark production 
based on a new perturbative resummation of gluon radiative corrections to the 
basic QCD subprocesses. We use Principal Value Resummation to calculate all 
relevant large threshold corrections. Advantages of this method include its in¬ 
dependence from arbitrary infrared cutoffs and specification of the perturbative 
regime of applicability. For pp collisions at center-of-mass energy yT = 1.8TeV 
and a top mass of 175 GeV, we compute cr{tt) = 5.52+q' 45 Pb- 



The quest for the top quark t reached fruition recently with the publication 
of results by two collaborations of experimenters studying tt pair production in 
proton-antiproton collisions at the Fer mil ab Tevatron 1 * 1 : p + p—>t + t + A\ The 
large mass of the quark and the possibility that its observed cross section exceeds 
theoretical expectations have, in turn, stimulated considerable theoretical activity. 
Among theoretical questions deserving attention is the quantitative reliability of 
cross sections based on the main production mechanism in perturbative quantum 
chromodynamics (pQCD), namely tt pair creation. 

At lowest order (tree-level), the two partonic subprocesses are quark-antiquark 
annihilation: 

q + q —> t + t (1) 

and gluon-gluon fusion: 

9 + 9 -> t + t . ( 2 ) 

These subprocesses are of order a 2 S in the strong coupling strength. The calculated 
top quark cross sections depend on these subprocess cross sections and on the 
parton density distributions that specify the probability densities of the quarks, 
antiquarks, and gluons of the incident p and p. 

Both order a ^ [21 and the next-to-leading order a jj [31 W contributions have been 
investigated thoroughly 151 . One observation is that the size of the 0(o? s ) terms 
in the partonic cross sections is particularly large near the tt production thresh¬ 
old, raising questions about the reliability of perturbation theory. This region 
of phase space is important for top quark production at the Tevatron owing to 
the large mass of the top quark. Virtual and bremsstrahlung corrections to the 
tree-level production channels, Eqs. (1) and (2), dominate the next-to-leading con¬ 
tributions in the near threshold region. The large threshold corrections can be iden¬ 
tified with numerically large logarithmic terms attributable to initial-state gluon 
bremsstrahlung, accompanying uncompensated mass-singularities removed after 
mass-factorization. After convolution of the partonic cross section with the parton 
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flux, the 0{a\) corrections to the physical cross section yield less dramatic en¬ 
hancements, of the order of 25% at top mass m = 175 GeV. These enhancements 
motivate a complete study of the large logarithmic corrections at the partonic 
level, to all orders in perturbation theory. In this Letter we present a resumma¬ 
tion of soft gluon corrections to the tt cross section employing the Principal Value 
Resummation (PVR) technique 161 . 

As in other hard-scattering processes, where large logarithmic threshold con- 
tributions are present , resummation of these corrections to all orders in a s is 
important both for theoretical understanding of the perturbative process and for 
numerical control of the resulting predictions. One method for resummation has 
been implemented previously for tt pair production W . In most resummation meth¬ 
ods 191 , including [8], the threshold corrections are exponentiated into a function of 
the QCD running coupling constant, a s , evaluated at a variable momentum scale, 
and formally integrated throughout the proper phase space. The variable scale is 
a measure of the momenta of the gluons emitted from the initial partons, and the 
associated integration represents the inclusiveness of the measurement. 

An inherent uncertainty and limitation of the resummation method of [8] is 
its dependence on an undetermined infrared (IR) cutoff. This cutoff enters when 
the exponentiated resummation function is integrated over phase space. In the 
formulation of [8], it cuts off the Landau pole of the QCD running coupling constant 
and the corresponding threshold region. Since the function is exponentiated, the 
dependence of the resummed cross section on this cutoff is important numerically. 
Further, dependence of the resummation exponent on this additional undetermined 
mass scale spoils explicit renormalization-group invariance properties for any region 
of the radiation phase space. Such properties can be implemented for a fixed value 
of the IR cutoff, but a simultaneous variation of the cutoff and the renormalization 
scale would again produce large sensitivities. 

The main advantage of PVR is that it does not depend on arbitrary IR cut¬ 
offs, as all Landau-pole singularities are by-passed by a Cauchy principal-value 
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prescription. We emphasize three strengths: 

(i) PVR reproduces correctly all the asymptotic properties and renormalon 
structure of pQCD. Due to the absence of extra undetermined scales, it allows for 
an evaluation of the perturbative regime, i.e., the region of the radiation phase 
space where perturbation theory should be valid. 

(ii) The approach implements universality of resummation automatically since 
it is the same mechanism for all relevant processes, independent of undetermined 
parameters. This mechanism produces a different perturbative regime for different 
production channels, but the mechanism itself is common to such channels. 

(iii) PVR been tested successfully in ll production [101 where its definition of 
resummed pQCD agrees fairly well with experiment, even at relatively low values 
of the hard scale (in the region of 5-10 GeV). 

In the remainder of this Letter, we present our results for the physical inclusive 
total cross section for ff-production in PVR for a top-mass range m £ {150,250} 
GeV, including a discussion of the remaining theoretical uncertainties. We also 
show some results for the behavior of partonic cross sections at m = 175 GeV. 
We present our predictions in the MS factorization scheme, in which the q , q 
and g densities and the next-to-leading order partonic cross sections are defined 
unambiguously. We postpone presentation of calculational details, more extensive 
comparisons with [8], and corresponding analyses of the differences to our longer 

[u] 

companion paper 

It has been observed that the functional form of the leading threshold correc¬ 
tions, order-by-order in perturbation theory, appears to be universal, i.e., indepen¬ 
dent of the hard-scattering process, except for differences in multiplicative color 
factors and process-specific kinematics. An example is the “universality” between 
dilepton- ( ll ) and ff-production, emphasized in [8]. We understand this universality 
to arise from the fact that mass-singularities in the former process are due purely 
to bremsstrahlung from the initial state, given that the ll pair is produced through 
an electroweak vector boson of large invariant mass, an “external potential” to the 
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QCD sector, whereas, in the latter the final quarks are (highly) massive, and gluon 
bremsstrahlung from (or interference with) that sector does not produce leading 
mass singularities and associated large logarithms. In tt production near threshold, 
the heavy quark final-state sector “factorizes” and behaves as an external heavy 
object, leaving the light initial states to produce the leading large QCD corrections. 
For the purposes of this work we accept this universality, assuming it valid to all 
orders in a s , as in [8]. Furthermore, we have chosen the order-by-order logarithmic 
structure and running coupling constant identical to those of [8]. 

The resummed partonic cross sections, including all large threshold corrections 
according to PVR, can be written as 

l 

a ij V (v, m2 )= J dz 1 + Hij(z,a) <r”(i/,m 2 ,z). (3) 

l-4(l+7/)+VT+i? 

In Eq. (3), 



crC(7/,m 2 ,z) = d{<r^\rj, m 2 , z))/dz, and is the tree-level partonic cross section 
expressed in terms of inelastic kinematic variables [8] to account for the e mi tted 
radiation* The variable 7/ = s/Am 2 — 1 is the distance from the partonic production 
threshold, a = a s (m)/7r, and ij £ {qq,gg} denotes the initial parton channel. The 
integration in Eq. (3) is over the phase space of the radiated gluons, parametrized 
through the dimensionless variable z. In ll production z = Q 2 / s, where Q is the 
invariant dilepton mass, and z = 1 corresponds to zero gluon momentum. 

Equation (4) is the main content of PVR. The leading large threshold correc¬ 
tions are contained in the exponent E{j(x,a), which is a calculable polynomial 

-k In other words, substituting 7i — > 0 in Eq. (3) and using the kinematic constraint 
crfj(r],m 2 , z = 1 — 4(1 + 77 ) + 4-y/l + v) = 0; we obtain the tree-level partonic cross sec¬ 
tion, < 7 ?-(? 7 , m 2 ) = < 7 ?-(? 7 , m 2 , z = 1 ). 
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in x. {Qj(x,a)} are calculable functions produced by the analytical inversion of 

the Mellin transform from moment space (where the argument of E{j corresponds 

to x —> Inn) to the physically relevant momentum space expressed in Eq. (4). 

These functions are produced by the resummation and are expressed in terms of 

successive derivatives of E^: Pk(x,a) = d k E(x,a)/k\d k x, k = 1, 2,.... Given that 

E contains at most one more power of x than of ct, Pf. contains at least k — 1 

fewer powers of x than of ct, and Qj contributes the resummed version of terms 

+ 

containing j fewer powers of x than of a in the integrand of Eq. (4). + For example 
Q o contains all powers PJ 71 , Q i all powers P 2 P™, etc. 

The main object of the resummation, Eij, was studied extensively in the con¬ 
text of ll production [10]. Embodying the universality discussed earlier, its func¬ 
tional form for tt production is identical to that of ll production, except for the 
identification of the two separate channels, denoted by the subscript ij. However, 
there are significant differences between the two processes which affect the way we 
implement PVR practically for the process at hand. These can be described as 
follows: 

First, only the leading threshold corrections are universal. The physical reason 
behind this is the existence of hadronic final states, absent in ll production, that 
produce interference effects with the initial hadronic states, a statement that can 
be deduced from the 0(a 3 ) calculation as well. Therefore, from all structures {Qj} 
in Eq. (4), the very leading one should be considered universal. This is the linear 
term in Pi, which turns out to be Pi itself. Hence, Eq. (4) can be integrated 
explicitly, and Eq. (3) may be written as 

l 

<7ij V (v, m2 )= j dzG Eii ^^' a ^a‘\ j {if,m 2 ,z). (5) 

l-4(l+7;)+VT+ii 

f For simplicity, we drop the channel indices. 

f These contributions are only approximately equal to the corresponding structures obtained 
upon a finite-order expansion of the resummation formula, since such expansion involves 
truncation of a product of series. 
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The second issue has to do with the perturbative regime in momentum space. A 
basic property of E{j (In n , ct) in moment space is that it has a perturbative represen¬ 
tation in the region 2b20tlnn < 1, where 62 is the first coefficient of the QCD beta 
function. In the region 1 < 262 ol Inn < 00 , it has a non-perturbative representation 
which, in turn, is an exponential suppression. In fact, lim^-nx, Eij(\nn, a) = — 00 . 
This latter “higher-twist” region is of negligible importance, especially for the large 
scale m of the top quark production process. These considerations apply in mo¬ 
ment space and are clearly independent of channel-dependent color factors, relying 
purely on the properties of the running coupling constant and PVR itself. To 
characterize a region in moment space as “higher-twist”, however, one must first 
convert to momentum space through inversion of the Mellin transform, Eq. (4). 
Specification of the boundary, mentioned already in ref. [ 6 ], is realized by the con¬ 
straint that all {Qj} , j > 1 be small compared to Qq which provides the leading 
integrand in Eq. (4), according to the previous power counting. This constraint 
can be shown quite generally to correspond to 

This perturbative regime in momentum space corresponds nicely with the perturba¬ 
tive regime in moment space, since Pi contains at most equal powers of logarithms 
and a. This determination of the perturbative regime is also less academic, in the 
sense that it realistically includes all relevant orders of perturbation theory, all rele¬ 
vant constants generated by the running of the coupling and, most importantly, the 
color factors that differentiate between various production channels. Color factors 
are not included in more abstract discussions in moment space but are physically 
extremely important. 

One could actually apply Eq. (5) all the way to z = 1 by using contour in¬ 
tegration as suggested by PVR, beyond the perturbative regime of Eq. ( 6 ), but 
one would then be using a model for non-perturbative effects, the one suggested 
by PVR, far beyond the knowledge justified by perturbation theory. In this work 
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we confine our attention to the perturbative regime, and hence our cross section is 
evaluated accordingly, namely 


a ij Vpei> (v,m. 2 ) = J dze Eii ^^-^a , ij {ri,m 2 ,z) 

l-4(l+7;)+VT+ii 


The upper limit is calculated through 


Pi In 


1 — Z r , 


a 1=1 


(7) 


( 8 ) 


Equations (7), ( 8 ) are the basic formulas of our approach. Based on these, 
notice that our perturbative approach can probe the threshold region down to 
77 > (1 — z max )l 2. It turns out our final result does not rely much upon the PVR 
method to by-pass IR renormalons and associated problems, precisely because it 
is restricted to the perturbative regime. In that sense, the presence of arbitrary 
IR cutoffs in previous resummations is superfluous, as all necessary information 
about IR sensitivity (i.e., the perturbative regime) can be obtained by examining 
the perturbative asymptotic properties of the resummation functions. Further 
analysis of these issues, especially the perturbative representation of Eq. (4) as 
an asymptotic series and related numerical issues relevant to the color factors 
characterizing each specific production channel, will be presented in our longer 
paperfll]. 

We provide our final resummed cross sections for each production channel by 
taking into account the complete next-to-leading calculation, through the 

improved prediction 


final/ 2\ PVpeit / 2\ , (0+1) x 2\ (0+1) /■ 2\ 

*ij (V,m) = *ij + (V,™ ) - a},- J (il,m ) 






(9) 


PV 


The last term in Eq. (9) is the next-to-leading order partonic cross section included 

PV 

in the resummation, a-- pert . 

> ij 
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In Fig.l we present our prediction for the partonic cross sections in the qq 
and gg channels, in the MS scheme, at m = 175 GeV. We also show the tree-level 
(cr( 0 )) and next-to-leading (cr( 0 + 1 )) counterparts. Notice that the three curves differ 
substantially in the partonic threshold region g < 1, with the final resummed curve 
exceeding the other two. Below g ~ 7 X 10 -3 in the qq channel and ~ 5 X 10 -2 in 
the gg channel, our final perturbative cross sections become identical to the next- 
to-leading order cross sections. This is a consequence of our decision to perform 
the resummation in the perturbative domain. This domain is smaller in the gg 
case because the color factors are larger in this channel. Above g ~ 1, we note 
that our resummed cross sections are essentially identical to the next-to-leading 
order cross sections, as is to be expected since the near-threshold enhancements 
that concern us in this paper are not relevant at large g. 

In Fig.2 we plot the physical cross section <r final , along with the tree-level and 
the next-to-leading order counterparts, for the qq and gg channels for a wide range 
of top mass. These are obtained after the partonic cross sections shown in Fig.l 
are convoluted with parton densities and integrals are performed over g. We use 
the most recent CTEQ3M parton distributions . Here, the usual factorization 
scale p has been chosen equal to the mass of the top. We observe that perturbative 
resummation enhances the cross section by about 10%. 

In Fig.3 we plot the same cross section as a function of the factorization 
(“hard”) scale p in a range /i/m £ {0.5,2}. Notice the mild dependence, as 
well as the interesting shape that peaks around the value 1. We consider the vari¬ 
ation of the physical cross section the range p/m £ {0.5, 2} a good measure of the 
theoretical perturbative uncertainty. 

Our prediction for the inclusive total ff-production cross section at the Teva- 
tron, using CTEQ3M parton densities, is 

i("> = U5 GeV) = 5,52±°;°J pb . (10) 

The “central” value of 5.52 pb is obtained with p/m = 1, and the upper and lower 
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limits of the uncertainty band correspond to the maximum and minimum values 
of the cross section in the range /i/m £ {0.5,2}. The cross section is insensitive 
to the choice of parton densities. Repeating the same analysis with the MRS(A') 
densities 1131 , we obtain 


= U5 GeV) = 5.321“; pb . (11) 


In Tables 1 (2) we summarize the main features of the physical cross section, 
in pb, for top mass m = 175 GeV for three choices of g and CTEQ3M (MRS(A')) 
parton densities. The first line in each table shows the values of the two-loop ol s / , k 
for the corresponding set; there is very little difference in the values. The MRS(A') 
cross sections are a bit smaller than the CTEQ3M values. For the qq channel this 
difference is less than 4% and is about 5% for the gg channel. 


m = 175 GeV 

fji = m/2 

g = m 

g = 2m 

Ct s /7T 

0.03584 

0.03260 

0.02991 

qq 

5.02 

3.69 

2.78 

qq 

4.13 

4.43 

4.21 

—final 

a qq 

4.70 

4.87 

4.55 

<T (0) 

a 99 

0.55 

0.34 

0.22 

a (o+A 

a 99 

0.72 

0.63 

0.50 

_.final 

*99 

0.75 

0.65 

0.52 

"final 

5.45 

5.52 

5.07 


Table 1. 
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m = 175 GeV 

p = m/2 

p = m 

p = 2m 

Ct s /7T 

0.03596 

0.03269 

0.02998 

a® 

_ <7 <7_ 

4.88 

3.57 

2.68 


4.00 

4.28 

4.08 

—final 

a QQ 

4.56 

4.71 

4.41 

n (0) 

a 99 

0.51 

0.32 

0.21 

n ( ° +1) 

a 99 

0.67 

0.59 

0.49 

—final 

*99 

0.70 

0.61 

0.50 

-tt 

"final 

5.26 

5.32 

4.91 


Table 2. 

The value in Eq. (10) is larger than that of [8], but within uncertainties, and 
it agrees with the present CDF and DO measurements [1], In Fig.4 we show the 
top-mass dependence of the physical cross section for pp —> ( tt)X. 

The band of perturbative uncertainty quoted in Eq. (10) is relatively narrow. 
On the other hand, we noted in discussing Fig.l that there is a reasonable range of 
7/ near threshold in which perturbative resummation does not apply. Perturbation 
theory is not justified in this region. Correspondingly, further strong interaction 
enhancements of the tt cross section may arise from physics in this region. We 
know of no reliable way to estimate the size of such non-perturbative effects and, 
therefore, cannot include such uncertainties in the estimates of the perturbative 
uncertainty of Eqs. (10) and (11). 

In our longer companion paper we will present mathematical and physical 
details, more extensive comparisons of our approach to existing results, and a 
discussion of remaining theoretical uncertainties. 
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FIGURE CAPTIONS 

1. The total parton-parton cross sections as a function of 7/ in the MS scheme 
at m = 175 GeV for the subprocesses (a) qq —> ttX and (b) gg —> ttX. Plotted are 
the lowest order Born cross section (dotted line), the next-to-leading-order cross 
section (dashed line), and the cross section obtained after resummation of soft 
gluons by the principal value resummation method (solid line). The QCD scale 
fi = m. 

2. The calculated cross sections for tt production as a function of top mass m 
in pp collisions at s/s = 1.8 TeV for the qq (upper bunch) and gg (lower bunch) 
subprocesses. Each bunch consists of three curves: the Born cross section (dotted 
line), the next-to-leading-order cross section (dashed line), and the cross section 
obtained after soft gluon resummation (solid line). 

3. Plot showing the calculated dependence of the final resummed cross section 
on (p/m) for tt production at m = 175 GeV and s/s = 1.8 TeV. Shown also is the 
next-to-leading order result (dashed curve). 

4. Physical cross section for pp —> ( tt)X at s/s = 1.8 TeV as a function of top 
mass. Data from the CDF and DO[1] collaborations are shown. Our calculated cross 
section here is obtained from adding the final resummed qq and gg contributions. 
Shown are entries for p/m = 0.5 (dashed), 1 (solid) and 2 (dotted). 
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cr (rj) (fib, gg, MSbar) 
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cr(m) (pb, i j , MSbar) 



m (GeV) 

Figure 2 
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/a/ m, m= 1 75 GeV 

Figure 3 
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